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ABSTRACT 



Velocity fluctuations over ocean waves 
are examined using joint probability densi 
tional mean function (JPDF-CMF) analyses, 
pared with theoretical predictions from wi 
models which provide for interaction betwe 
duced motion and turbulence in the airflow 
predictions from potential flow theory. 

Phase and amplitude results for wave-r 
tions are examined in detail and consisten 
potential flow are noted. These departure 
ences in between u and w which are not in 
tudes of u and w which vary throughout the 
phase shifts of *u and w which differ from 
by potential flow theory. 

Conclusions for the role of turbulence 
tions are based on JPDF-CMF and phase ampl 



and wave heights 
ty function-condi- 
Results are com- 
nd-wave coupling 
en the wave in- 
, and also with 

elated fluctua- 
t departures from 
s are phase differ- 
quadrature, ampli- 
per iods , and 
those predicted 

in the interac- 
itude analyses. 
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INTRODUCTION 



More detailed studies of the near atmospheric boundary 
layer structure and energy exchanges are necessary if solu- 
tions to several applied geophysical problems are to be 
solved. Some of these applied geophysical problems include: 
the forecasting of surface wave conditions which require 
wind speeds at specified layers in the surface layer for 
semi-empirical formulae, and specifications of kinetic 
energy losses and thermal energy gains in the surface layer 
that are necessary for extended period weather predictions. 
Previous methods for obtaining boundary fluxes and deter 
mining the structure of the air layer near the ocean surface 
consisted primarily of applying empirical relationships, 
which have been found to be valid over land. These previous 
efforts generally made use of the assumption that interac- 
tions associated with wind-wave coupling are negligible. 

'An example of the need for a better understanding of the 
atmospheric boundary layer structure was the discrepancies 
in the interpretations of results presented at the Confer- 
ence on the Interaction of the Sea and Atmosphere at Fort 
Lauderdale, Florida in December 1971. Various results were 
marked by conflicting interpretations when examined with 
strictly empirical approaches. Most evident in these were 
the various representations and values obtained for the drag 
coefficients as well as differences from in ter comparisons 
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between stress estimates obtained from different methods. 
The basic assumptions in obtaining these results were the 
form of the wind profile and neglect of significant depend- 
ence on the waves. General descriptions of the influence 
of wind-wave coupling within the framework of theoretical 
models combined with observational data should yield re- 
sults which delineate reasons for these discrepancies. 

It is appropriate to perform detailed analysis of velo- 
city fluctuations over natural waves at this time not only 
because of the need to consider the dynamics of wind-wave 
coupling to explain discrepancies among results but also 
because of the present stages in both theoretical and obser 
vational efforts in this area. Recent advances in theoreti 
cal and observational studies have brought both to the 
point where the results should and can be related. 

Theoretical studies on wind-wave coupling have pro- 
gressed from the initial quas i-laminar descriptions (see 
Miles 1957) wherein the wave-induced motion interacts only 
with the mean wind to those wherein interaction occurs with 
the turbulence in the airflow (Yefimov 1970). Solutions 
from theoretical models which allow for turbulence in the 
airflow have revealed several features recently observed in 
results obtained over natural waves. These features in- 
clude amplitude and phase variation of the wave-related 
fluctuations with height and contributions by the wave- 
induced motion to the cospectra of u and w. Similarities 
in these results suggest that theoretical models could be 
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refined and the observational results interpreted better if 
the two types of results were inter compared • 

In this study, velocity fluctuations over natural waves 
will be examined in order to show that the airflow does not 
necessarily follow the potential flow definition inherent in 
early wind-wave coupling models* For this, recent theories 
and discussions of wind-wave coupling will also be con- 
sidered. These include those by Yefimov (1970), Manton 
(1971), Davis (1972), and Davidson (1970). 

The data examined represent simultaneous velocity mea- 
surements at different levels over ocean waves on the Atlan- 
tic Ocean. The measurements were made from the Scripps 
Institution of Oceanography Floating Instrument Platform 
(FLIP) which appears in Figure 1. The measurements were 
made during the Barbados Oceanographic and Meteorological 
Experiment (BOMEX) . The BOMEX experiment was a joint inves- 
tigation of several agencies and universities for the pur- 
pose of establishing the energy budget of an atmospheric 
volume 500 km on a side and about 5 km high. The volume was 
located in the Trade Wind region of the Atlantic Ocean as 
shown in Figure 2, The data represent measurements near the 
center of the volume, in the vicinity of 16°N latitude and 
57°W longitude. 

In this study, joint probability analysis on velocity 
and wave are the primary source for information on inter- 
actions that could have taken place due to wind-wave 
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Figure 1. Side and Top View of FLIP. 
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Figure 2. BOMEX and FLIP Location During 19-28 May 1969. 
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coupling. From the analyses, interpretations will be made 
on phase-amplitude relations, non-linear properties and 
determination of the type of turbulence present. It is sug- 
gested that these statistical procedures will find wide- 
spread acceptance in the study of micro-scale meteorological 
phenomena . 

Data examined in this study were also analyzed by spec- 
tral analyses in a parallel investigation by Bingham (1972) . 
Some of his results will be referenced during the discus- 
sions. Spectral analysis, however, is limited when applied 
to studies involving phase-amplitude relationships between 
wind and wave. 
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II. ASPECTS OF THEORETICAL AND OBSERVATIONAL RESULTS 



A primary motivation for the design of this investiga- 
tion was the close agreement between recent observational 
results and numerical solutions from recent wind-wave coupl- 
ing models by Yefimov (1970) and Reynolds (1968). In for- 
mulating their models, Yefimov and Reynolds allowed for 
interaction between the wave-induced motion and the turbu- 
lence in overlying shear flow. Although most theoretical 
formulations of the linear theory for wind-wave coupling 
have several common assumptions; such as, two dimensional 
flow and a specified mean wind dependent only on height, 
and common simplifications; such as a single component wave 
field and suitable linearized forms for the defining equa- 
tions, a dis ting*uishing feature of recent formulations has 
been the role prescribed for turbulence in the airflow. In 
contrast. Miles (1957), in formulating an early quasi- 
laminar model, neglected turbulence except for its role in 
determining an assumed logarithmic mean profile. 

Davis (1972) has shown that anistropic stresses, which 
are most likely to be turbulent Reynolds stresses, are 
essential to the process of energy flow to the wave. These 
stresses, however, have been difficult to isolate in the 
various experimental tests of the models. Careful measure- 
ments and improved analyses of the fluctuating velocity 
components and the turbulent stresses are apparently what 
are required . 
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A, DISCUSSION OF WIND-WAVE COUPLING THEORY PREDICTIONS 

Miles* theory is a quasi-laminar wave generation mechan- 
ism in which the wave-related motion is assumed to follow 
the potential flow theory. The key to Miles* formulation 
is the existence of the critical level which corresponds to 
the level where the wind speed is equal to the wave speed. 

A schematic of Miles* assumed airflow is shown in Figure 3. 
The wave in this figure is stationary. Below the critical 
level (dotted line) , the air appears to move to the left 
since it moves at a slower speed than the waves. Above the 
critical level, the wind moves to the right. Miles* theory 
is based on the premise that waves only modulate the mean 
wind flow over them and that turbulence is not influenced. 
Reynolds stresses are then dependent only on the height 
above the mean wave surface.* Phase and amplitude features 
of the two components of velocity, and the stress uw, with 
respect to the wave, as predicted by potential flow theory 
are illustrated in Figure 4. 

Recently, Yefimov (1970) and others have allowed for 
turbulence in the airflow, which requires taking into 
account nonlinear terms in the equations of motion. A way 
to approach this has been to define the total vertical 
transfer of horizontal momentum as, (Yefimov 1970), 

-u'w' = ^ + Y(tT 

m o z 3 z o z 

where u* and w* represent total velocity fluctuations, u and 

w are wave-related fluctuations, u is the mean wind, k = 

m 
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Figure 3. 



Streamlines of Sinusoidal Perturbations 
in a- Shear Flow Near the Critical Level. 
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Figure 4. Phase-Amplitude Relations Predicted by Potential 
Flow Theory for the Two Components of Velocity 
(u and w) and the Stress uw, with Respect to the 
Wave . 
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ku^z and y is an exchange coefficient related to the wave- 
induced motion. 

Several features in Yefimov’s solutions are due solely 
to "the prescribed interaction between the wave-induced 
motion and turbulence. As such, they occur independent of 
the height or even existence of a critical level. In com- 
parison, all nonpotential features of the airflow in initial 
quasi-laminar formulations were intimately connected to the 
dynamics at the ’critical level’. These features were the 
variations with height of the phase difference between the 
velocity components and the wave. Phillips (1966) and 
Kinsman (1965) describe the predicted phase shift for the 
vertical velocity component as being from over the downwind 
node toward the crest as the ’critical level’ is approached 
from the surface. Stewart (’1969) and Reynolds (1968) were 
able to verify results from a wind-wave tunnel experiment 
with phase predictions from a nonlinear formulation for 
wind-wave coupling and their work will be referenced for 
this aspect. 

Consideration of results with respect to theoretical 
approaches for the consequences of wind-wave coupling will 
be confined to prescribed or predicted characteristics of 
the velocity field in these models. These characteristics 
include ; 

a. Variation of the amplitude of wave-related fluctua- 
tions with height; 
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b. Contributions (negative or positive) to the 
momentum flux; and 

c. Variation of the phase difference between wave- 
related velocity .fluctuations and the wave with 
height . 



B. STATISTICAL CONSIDERATIONS AND PROCEDURES 

Both joint probability density function-conditional mean 
function analyses (JPDF-CMF) and spectral analyses will be 
used to identify properties of fluctuation velocities asso- 
ciated with the waves’ presence but only the JPDF-CMF analy- 
ses will be described in some detail. 

The following diagram outlines procedures^ used in the 
JPDF-CMF analyses. The procedures will be discussed in the 
ensuing paragraphs . 



Normalization and Computation 
of Preliminary Statistics of 

the Variable 
_ 



Band Pass Filtering f 




1 . Preliminary Statistics Computation and Normalization 
The procedures were applied to normalized values of 
the original data. Normalization was accomplished, by 



All statistics were compiled and plotted at the W. R. 
Church Computer Facility of the Naval Postgraduate 
School utilizing the IBM computer and a CALCOMP plotting 
device . 
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dividing the deviation from the mean by the standard devia- 
tion; for example (x - J O * 

2 . Band Pass Filtering 

A band pass filter was applied in order to focus 
the analyses on fluctuations near the frequency band corres- 
ponding to predominant surface waves. Specifically, 
application of the band pass filter allowed analyses of 
wave-induced fluctuations to be clear of any high or low 
frequency contamination, A numerical inverse transform 
filter was applied and was designed on the basis of pre- 
vious spectral results. The response curve for the selected 
band pass filter appears in Figure 5, where low and high 
frequency cut-offs are observed to be approximately ,1 and 
,8 Hz respectively. 




Figure 5, Numerical Inverse Transform Filter, 
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3. 



JPDF-CMF Computations 



JPDF-CMF procedures adopted for this study have been 
fully described by Holland (1972). In general, a JPDF re- 
presents the probability of a joint occurrence for a pair 
of variables and the CMF represents the average value of a 
third variable in a joint class interval for all observa- 
tion times when JPDF variables had values within the joint 
interval. Figure 6 shows a two-dimensional array represent- 
ing the trivariate statistical relationship described by 
the JPDF and CMF. From this figure, the probability of 
joint occurrences in two variables can be compared with a 
corresponding conditional mean of a third variable. 

4 . Phase-Amplitude. Theory and Computation 

Available theoretical (Yefimov 1970, Manton 1971) 
and observational (Davidson ' 1970 , Holland 1972, Davidson 
and Frank 1973) results indicate that one of the important 
features of interest in this type of study is the statisti- 
cal dependence of fluctuations in the airflow on the phase 
of the wave. Such relations were determined by Davidson 
and Frank (1973) using procedures which are used in this 
study. It is in regard to phase-amplitude information that 
greater latitude for interpretations can be made in JPDF- 
CMF analyses than those obtained by spectral analyses. For 
example, JPDF-CMF results should reveal nonlinear features 
since the information is not a result of linear expansions 
of the fluctuations. Furthermore, the approach can be 
applied to comparisons between wave-amplitude classes. In 
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Figure 6. Two-Dimensional JPDF-CMF Array. 
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this study the features are examined with respec 
large, and combined or total waves. 

Phase information was determined using a 
cribed by Holland (1972) which is based on JPDF- 
procedures. Phase-amplitude information is poss 
such procedures because polar coordinates which 
the phase angle and amplitude of a variable such 
wave height (r|) can be determined from the JPDF 
the variable, f| , and its time derivative, f| . As 
of this, u (fj , ri) would show the dependence of th 
tal component of wind velocity on the phase and 
of the waves. As will be seen in later discussi 
phase is defined in terms of octants measured co 
wise in the fi , X] coordinate system (JPDF array) 
positive T*| axis.* JPDF-CMF arrays which contain 
amplitude information for u,v,w,uv,uw, and vw wi 
to the waves will appear in the discussions of a 
section. The phase-amplitude results are easily 
by octants for each period analyzed. 

In determining phase-amplitude informati 

meter, T^, called the 'statistical eddy period' 

interest. t is defined by the ratio 
o 

T = 2tt ~ 

° d 

where a and 6 are standard deviations of the var 
its derivative. If the statistical eddy periods 
differ greatly between the independent and refer 
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method des- 
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variable, phase-amplitude informat 
variable has a greater possibility 
mean of numerous random fluctuatio 
statistical eddy periods agree as 
With respect to possible i 
on the results, analyses by Rudnic 
that observed fluctuations in u an 
than could be induced by FLIP’s mo 
motion will not enter into the dis 



ion for the reference 

of only representing the 

ns. In general, the 

can be seen in Table I. 

nfluence by FLIP’s motion 

k (1967) have indicated 
2 

d w are 10 times larger 
tion. Therefore, FLIP’s 
cussion of the results. 
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Table I. Observed and Predicted Velocity Components, and 
Statistical Eddy Periods for Each Variable. 

Period U-Component (M/Sec) W-Component (M/ Sec) Statistical Eddy Period(t 

Subset Observed Predicted^ Observed Predicted^ n u v w ■ uv uw 
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DESCRIPTION OF MEASUREMENTS 



The data available for this study consisted of simul- 
taneous recordings of wind velocity components (with hot- 
film anemometers) , air temperature fluctuations and wave 
heights. The levels .of measurement ranged from 2 to 8 
meters above mean water level and the sensors and recording 
systems were capable of measuring velocity and temperature 
fluctuations of frequencies up to 50 Hz, There were about 
54 hours of recordings in 40 separate observation periods 
to select from. On the basis of meteorological and wave 
conditions and on previous spectral analyses, four periods 
were chosen for the analyses. 

In the following paragraphs specific details on the 
measurements and general conditions for every period will 
be described, 

A, FLOATING INSTRUMENT PLATFORM (FLIP) 

A discussion of FLIP has been provided by Bronson and 
Glosten (1968), The Floating Instrument Platform is a 
super-stable, open-sea, free-floating platform, which is 
towed to a site while in a horizontal position and then 
partly flooded to cause the structure to flip to a vertical 
position for use as a research platform. Figure 1 included 
a side and top view of FLIP, while Figure 7 illustrates the 
working section of FLIP which remains above the water line 
when FLIP is in the vertical position. 
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Figure 7. Sensor Locations on FLIP. 
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In general, FLIP performed in a satisfactory manner 
throughout its use in BOMEX, Wind tunnel tests revealed 
that a minimum distortion to the natural flow of air past 
the intended location of an outboard mounted vertical sen- 
sor supporting mast occurred when the deck of FLIP was 
positioned facing into the wind. 



B, EQUIPMENT AND PROCEDURES 



1 • Sensor Mounting Arrangement 



Figure 7 shows the 


locations of 


sen 


sors for 


turbu- 


lence, profile and 


v;ave me 


asurements on 


the 


ver t i cal 


mas t 


attached to FLIP. 


Sensors 


were mounted 


at 


the 2 , 3 , 


4 , 6 , 


and 8 meter levels 


and all 


sensors were 


more than 1 


meter 



windward from the vertical mast. As can be seen in Figure 
7, the vertical mast, supported by a horizontal boom, is 
about 15 meters crosswind from FLIP's hull. The mast, was 
attached to a carriage that could move horizontally on the 
boom. In addition, the mast could be turned on a horizon- 
tal axis so that it could be brought alongside a working 
platform to replace and adjust sensors during the experi- 
ment , 



No provision for adjusting the 
of the mast, independent of FLIP, were 
made it almost impossible to maintain 
level position. Provisions were made, 
for this discrepancy. In general, it 
ment mounting arrangements have not si 



ver tica 
availab 
sensors 
however 
is felt 
gnif ican 



1 position 
le. This 
in a mean 
, to’ correct 
that instru- 
tly affected 
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the accuracy of the data values with regard to 
tations made. 

2 . Wind Component Measurement 

Wind components were measured with hot 

2 

stant temperature anemometer systems, with lin 
3 

Each probe had three quartz-coated mutually p 
sensors 0.15 millimeter in diameter and 2.0 mi 
long. Each sensor is a glass rod covered with 
film onto which a quartz coating has been sput 
Figure 8) . The sensors were electrically isol 
simultaneous measurements could be made. The 
array made it possible, with trigonometric rel 
to determine the instantaneous vector componen 
wind from simultaneous measurements. The freq 
of these sensors, appeared to be well over 50 H 
conditions under which they were used. 

The constant temperature anemometer sy 
velocity by sensing the heat loss from the hea 
to the moving air. If the air density is near 
and if the sensor is maintained at a constant 
ly high temperature, the heat loss is dependen 
entirely on air velocity. The sensor is part 
circuit . 
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Figure 8. Three-Sensor Hot Film Probe. A Millimeter Scale 
is Shown on the Left. 
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3, Wave Measurement 



Wave data were obtained from a resistance gauge 
made available by Dr. R. E. Davis, Scripps Institution of 
Oceanography. The probe consisted of a non-conducting 
tube, approximately 2.5 centimeters in diameter, with a 
conducting wire wrapped spirally around it. It was posi- 
tioned about 5 meters inboard from the vertical mast that 
held the velocity and temperature sensors. 
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IV. PRESENTATION OF RESULTS 



The results are from four different periods. The per- 
iods were in three separate days in May 1969: 19 May (1613 

-1735 GMT), 27 May (0311-0435 GMT), and 28 May (1756-1916 
and 2028-2150 GMT) . Each period was approximately 80 min- 
utes in duration. Hereafter, these periods will be referred 
to as periods 1, 2, 3, and 4. Each period consisted of 
several subsets with regard to JPDF-CMF and phase-amplitude 
analyses; periods 1 and 2 consisted of four subsets of 20 
minutes each, and periods 3 and 4 consisted of three sub- 
sets of 20 minutes each. Separate reference will be made 
to individual subsets within these periods when required. 

Airflow data for periods 1, 2, and 4 were from the 8 
meter level and those for period 3 were from the 3 meter 
level. Data considered were the u, v, and w velocity com- 
ponents and wave heights. 

General conditions during the four periods are sum- 
marized in Table II and Figure 9. General conditions 
during all periods could be described as being steady with 
respect to both waves and wind. Certainly the periods did 
not have active wave growth or rapid decay. These condi- 
tions are representative of conditions over most of the 
oceans when swell propagates into regions of relatively 
light winds. Therefore, these results are relevant to des- 
cribing the surface layer over the ocean. 
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Spectral results from these periods have been described 
by Bingham (1972), They show corresponding extrema in the 
velocity variance, covariance spectra and wave spectra and 
reflect a few characteristics of the predicted solutions, 
as will be noted in the results. UW cospectra had extrema 
which, in fact, represented either enhanced or decreased 
momentum transport when compared to that at frequencies not 
associated with wave spectra extrema. Potential flow does 
not predict such extrema in the uw cospectra because in 
theory, u and w are in quadrature. Such contributions, 
both negative and positive, are interpreted as being those 
predicted by the recent theoretical solutions. 

Physical processes can be realized from the arrays re- 
presenting probabilities of joint occurrences of fluctua- 
tions of various magnitudes in (u, w, uw , etc.) in 
conjunction with an average surface elevation during the 
joint occurrences. 

JPDF-CMF patterns for potential flow are presented in 
Figure 10. As indicated in Figure 4, potential flow 
theory predicts a maximum positive value of uw occurring 
in quadrant II and IV. The JPDF-CMF pattern in Figure 10, 
(Panel A), represents such a relationship. Similarly, u 
and w relationships which appear in Figure 4 are repre- 
sented in Panels B and F. Panel C and E will be discussed 
in the following paragraph. Panel D reflects the absence 
of any correlation between the CMF , v, and the JPDF (fi, n) • 
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TABLE II 



General Conditions for Periods Analyzed 



PERIOD 


MEASUREMENT 

LEVEL 

(METERS) 


WIND 

SPEED 

(M/SEC) 


PHASE 

SPEED 

(M/SEC) 


CRITICAL 

HEIGHT 

(METERS) 


(M/SEC) 


C 

u* 


1 


8 


9.7 


13.9 


> 


15 


.20 


69 


2 


8 


11.6 


16.1 


> 


15 


.28 


58 


3 


3 


10.5 


15.2 


> 


15 


.17 


88 


4 


8 


9.3 


14.7 


> 


15 


.24 


62 
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Panel C represents the probability of joint occurrences 
of fluctuations of various magnitudes of u and w with res- 
pect to the wave height. The following diagram interprets 
the four quadrants in Panel C with respect to momentum 
transfer. Panel E has, of course, the CMF , wave derivative, 
rotated 90 degrees clockwise from its orientation in Panel 
C. 



Upward transfer of 
deficit momentum 

-U 



+ W 

Upward transfer of 
excess momentum 

+U 



Downward transfer of 
deficit momentum 

- W 



Downward transfer of 
excess momentum 



The following sections include detailed discussions of 
features in JPDF-CMF patterns for each of the four periods 
commencing with period 1 and continuing through period 4. 

A summary will follow the presentation of results. 

A. PERIOD 1 

Significant features identifiable within JPDF-CMF re- 
sults from this period are evidence of mechanical turbulence 
associated with hydrostatic stability, decreased momentum 
transfer due to the waves’ presence, large deviations from 
potential flow and apparent changes in u and w phase rela- 
tions with respect to the wave due to height variations of 
the critical level. 



35 



mAMB 









m 






















o 


CO 














O 












CO 


4J 












dJ 


Cd 














e 












•H 














4J 
















CO 












n:J 


0) 












a 














cd 


1H 














cd 












CO 


* d 












4J 


a CO 














M -H 












W) 1-4 > 












*H 






• 








(U 


1 




CO 






• 




'd 




u 






w 




M I 




0) 






a 




Cd 1 




4J 






o 


o 


O 1 




Cl) 






•H 




1 










U 


CO 


Cd 1 










•H 








o 








CO 


CO 




rH 




h- 


a 




4J CO 








2 


o 


CD 


C 0) 




4J 




O 


o 


4J 


0) > 




Cd 








4J 


6 cd 








LU 

L-, 


o 


(U 


o ^ 




•X? 




> 


tH 


u 




0) 




< 






d I 




0) 




Q 






CO 




a 








c 


cd X 




CO 








cd 


Q> 






• 




c 




n - 




no 






c3 


CO 






d 


O 






rH 


d rH 




•H 








O 


Cd d 






a\ 






rO 


to O 






H 




•H 


e 


•H 




0) 








H 




bO .H 






CO 


O 




cd 


Cd 








•H 1 


• 








C 


Q) 




CO 


O 


4J 




oj 




H 


4J 


> 


0) 






H 


4H 




<C 






w 




o • 


bO 




d 




4J 


• • 


H 


•H 


• • 


Cd 




d 


U-l 




0) 


4H 


e 




OJ 


rH 


4J :$ 




rH 






e 


cd 


•H ^ 




Cd 


U 




0) 




CO > 


0) 




O 




M 






> 




P 




:3 




(1) d 


cd 








w 


Cl) 


> 


D5 


0) 


6 




otf 


a 


•H 1 




:5 


o 




o 


a 


d 


4H 


o 


u 




:si 


P| 


P • 


O 


iH 


Pm 



' — ' 



1H0I3H 






U 




Figure 10. JPDF-CMF Patterns for Potential Flow. 
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With regard to the evidence of mechanical turbulence, 
the JPDF (u,w) patterns can be interpreted for these fea- 
tures in the fluctuations, Holland (1968) proposed that if 
turbulence is primarily mechanical and is occurring under 
stable conditions then a wider distribution of both u and 
w should occur at positive u than at negative u and the 
main mode, at high probability should be shifted towards 
negative u. Distributions which fit such a description 
appear in the JPDF pattern in Figure 11, Panel C. If the 
hydrostatic conditions had been unstable, the opposite 
effect, a narrow distribution of u and w at positive u and 
the main mode shifted tov/ards positive u would have been 
seen . 

JPDF patterns below the 5 percentage level in Figure 11, 
(Panel C) are slightly elliptic, asymmetrical with a major 
axis 135 degrees counterclockwise from the positive u~axis. 
Statistically, a small probability of joint occurrences of 
fluctuations of various magnitude of u and w is observed 
in quadrants II and IV, indicating that, due to the waves’ 
presence, upward transfer of deficit momentum and downward 
transfer of excess momentum occurred. 

From cospectra in Figure 12, it is observed that upward 
transfer of momentum was evident near the wave spectra 
peak. Similarly, the center patterns in Figure 11 (Panel 
C) are nearly circular and correspond to decreased downward 
transfer as shown in the cospectra. Figure 12, Recall that 
in these interpretations the band-pass filter (see Figure 5) 
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Figure 11. JPDF-CMF Patterns for Period 1, Subset 1. 
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included frequencies of a greater range than the wave 
spectra peak band. 

With respect to phase relations, large deviations from 
potential flow can be noted by comparing u(fj, n) and w(n, 
ri) patterns of Figure 13, (Panels B and F) with comparable 
panels in Figure 10. In particular, a rotation of about 
45 degrees is seen in the CMF of u(ti, T]) . This is an 
obvious shift from that predicted by potential flow (Figure 
10, Panel B) . CMF patterns in w(f|, ri) are not as dominant 
as the patterns to be considered later for period 3, yet a 
rotation of about 45 degrees can be discerned. The phase- 
amplitude relationship in Figure 14 further reflects the 
above deviation from potential flow. The phase-amplitude 
representation shows w to be nearly in-phase with the u 
component. In contrast, potential flow theory predicts 
that u and w will be in quadrature as was shown in Figure 4. 

Karaki and Hsu (1968) observed phase results similar to 
those just described by comparing u and w perturbation over 
water waves in a laboratory experiment. These were ob- 
served, from their analysis, near the region of the 
critical level and were due to a shift in the phase of w 
relative to both the wave phase and u-phase. The relative 
u-w phase was observed to be near zero in the critical 
level. Stewart (1970) also described results where the 
phase shift of both u and w with height resulted in u and w 
being out of quadrature and nearly in-phase. 
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Figure 12. Spectral Results 1) Including Covariance 

Spectra and Wave Spectra. 
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Figure 13. JPDF-CMF Patterns for Period 1, Subset 3. 
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Figure lA . Phase^Amplitude Results (Components) for 
Period 1, Subset 2. 
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The relations just described can be compared to those 
occurring later in the period by comparing the results in 
Figure 14 to those in Figure 15. Results in Figure 15 cor 
respond to a time interval which was 42 minutes after the 
interval corresponding to the results in Figure 14. 

Considering Figure 15, the phase difference between u 
and w gradually increased as the period progressed until a 
shift of nearly 45 degrees was reached during the last sub 
set of the period. Figure 15. 

B. PERIOD 2 

Significant features observed in this period were: 
evidence of mechanical turbulence associated with hydro- 
static stability, apparent potential flow conditions and 
variations in the magnitude and direction of the wave-in- 
duced momentum transfer. Also changes appear to have 
occurred during the period in the physical make-up of the 
wave. This period, unlike any other period, had near 
potential flow conditions in both u and w. 

With respect to the role of mechanical turbulence, in 
Figure 16 (Panel C) ri(u,w) JPDF patterns are highly ellip- 
tic with asymmetrical rings orientated so that the major 
axis is 135 degrees from the positive u-axis . The highly 
elliptic pattern will not be observed in any of the other 
periods. Again the turbulence is .interpreted to be 
mechanical as in period 1 on the basis of the shift of the 
JPDF pattern. The corresponding cospectra in Figure 17 
show only slight enhanced momentum transfer. 





Figure 15. Phase-Amplitude Results (Components) for 
Period 1, Subset 4. 
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JPDF-CMF Patterns for Period 2, Subset A. 
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Spectral Results (Period 2) Including Co- 
variance Spectra and Wave Spectra. 
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Figure 17 . 




9 



♦ 






i 




With respect to similarities with potential flow pre- 
dictions, u(f|,n) w(f),n) results in Figure 16 (Panels 

B and F) and throughout the period reveal no significant 
shift of the CMF patterns from those predicted for poten- 
tial flow (Figure 10). Also, w(n,n) CMF patterns are weak- 
ly defined when compared to other periods. This indicated 
a smaller probability of joint occurrences in average 
fluctuations of w with f\ and n* 

With respect to the suggested possibility of a differ- 
ent wave shape, JPDF patterns in Figure 16 (Panels B, D, 
and F) have a marked difference from those in any other 
period, such as those in Figure 11. Where previous and 
later periods have nearly concentric patterns showing only 
slight deviations from those of a sinusoidal wave, this 
period has mor e ellip t ical patterns with a slightly in- 
creased gradient over the wave crests. The corresponding 
’statistical eddy period’ is also much lower during this 
period than the other periods (Table I). This could be 
indicative of a marked change in the physical make-up of 
the wave, waves having sharper crests and broader troughs. 

With respect to momentum transfer, uw , in relation to 
amplitude classes, the large wave results in Figures 18, 

19, 20, and 21 appear to shift from an upward transfer of 
deficit momentum to a downward transfer of deficit momentum 
as the period progresses. To a lesser degree, the shift 
is also present in the total wave. This is interesting in 
view of the results in Figure 4 which show the general wind 
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Figure 18. Phase-Amplitude Results (Stresses) for 
Period 2, Subset 1. 
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Figure 19. Phase-Amplitude Results (Stresses) for 
Period 2 , Subset 2. 



50 



I 




♦ 



« 





Figure 20. Phase-Amplitude Results (Stresses) for 
Period 2, Subset 3. 
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Figure 21. Phase-Amplitude Results (Stresses) for 
Period 2, Subset 4. 
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conditions for all periods. The wind has been increasing 
with a large increase during the last part of this period. 
This can possibly be attributed to what Yefimov (1970) 
called the complex phase as the wind speed approaches the 
phase speed of the predominant surface wave. 

The last feature noted for this period is that the wave- 
form of the wave related stress (uw) , Figures 18, 19, 20 
and 21, also varies considerably over the period. Consider- 
ing the large wave amplitude class, early in the period the 
largest amplitude exists in quadrant II but as the period 
progresses, the dominent quadrant shifts to quadrant I. 

With respect to the total waves, however, initially the 
wave-form is similar to the potential flow prediction but 
quadrant II becomes dominant later in the period. With res- 
pect to the third subset, the large wave amplitude class 
has the maximum value in quadrant III as does the total 
wave amplitude class. The final minutes of this period 
continue with the same relation found in subset three, that 
of maximum amplitudes in quadrant III. 

C. PERIOD 3 

This period differed from the other periods in that 
measurements were made at 3 meters. Significant features 
in period 3 appear to be reduction in excess momentum 
transfer at the frequency of the wave spectrum peak, dis- 
tinct deviations from potential flow in the phase of w, and 
uw phase-amplitudes which differ significantly from potential 
flow theory. 
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With respect to momentum transfer, JPDF patterns in 
Figure 22, (Panel C) have the appearance of concentric cir- 
cles. They do have a slight asymmetry in quadrant IV, re- 
presenting some downward transfer of excess momentum. 
Spectral results for this period appear in Figure 23. The 
uw cospectral results reflect excess momentum transfer at 
the frequency of the wave spectrum peak. Such a spectral 
result does not agree with the JPDF axis rotation in Panel 
C, which is approximately 315 degrees counterclockwise from 
the positive u-axis . 

The CMF patterns are now examined in order to determine 
the role of the waves in producing the (u,w) JPDF patterns 
in Figures 22 and 24 (Panel C) . For potential flow the 
CMF patterns for the waves would have maximum positive 
values associated with negative u and zero w (Figure 10, 
Panel C) . Maximum negative CMF values would be associated 
with positive u and zero w. A distinct deviation from the 
pattern is the clockwise rotation of both maxima and minima 
in the CMF patterns of approximately 35 degrees in Figures 
22 and 24 (Panel C) . 

With respect to phase relationships, the significance 
of the phase shift in both extrema can be evaluated further 
by examining the JPDF-CMF patterns for u(f|,ri) and w(f|,ri) in 
Panels B and F. The CMF pattern orientation predicted by 
potential flow describes that for u fairly close with only 
a slight rotation as shown in Figure 22, Panel B. In the 
case of w, however, where maxima/minima extremes are 
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Figure 22. JPDF-CMF Patterns for Period 3, Subset 1. 
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Figure 23. Spectral Results (Period 3) Including Co- 
variance Spectra and Wave Spectra. 
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Figure 24. JPDF*-CMF Patterns for Period 3, Subset 3. 
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supposed to be concentric around zero values of r\. There rs 
a distinct rotation backward toward positive and negative 
values of n- This also can be seen from the summarized 
phase-amplitude for this period, Figure 25. For all ampli- 
tude classes, u is also comparable to potential flow, but 
w is not. The latter is made evident by the minimum value 
of w occurring before the inflection point on the wave. 

This represents, perhaps, evidence of downward transfer of 
excess momentum over the phase of the wave due to 
mechanical turbulence. 

To further emphasize the distortions that appear in 
Figure 25 for u and w, the uw phase-amplitude results and 
uw(f|, n) JPDF-CMF results will be considered. In potential 
flow. Figure 4, maximum values of uw occur behind the crest 
and following the trough, whereas minimum values occur 
immediately preceding the crest and the trough. The phase 
relations appearing in Figure 26 are comparable to those 
predicted; however, the amplitudes vary over the different 
phases of the wave. For example, quadrant IV amplitude 
values are greater than the others. This is also apparent 
in the JPDF-CMF results. The latter result in some down- 
ward transfer of excess momentum. Quadrant III ampli- 
tudes are also greater than those in quadrant I which 
suggests that the occurrence of downward transfer of 
deficit momentum predominates as associated with hydro- 
static stability, and, therefore, consistent results are 
obtained in this period. 
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Figure 25. Phase-Amplitude Results (Components) for 
Period 3, Subset 2. 
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Figure 26, Phase-Amplitude Results (Stresses) for 
Period 3, Subset 1. 



60 






I 



I 



I 

< 



Considering the same features as above but for the JPDF- 
CMF uw(f|, ri) , this period has maximum values of uw, Figure 
22, (Panel A), occurring in quadrant IV but rotated slightly 
towards the trough position. This rotation, although not 
as apparent in the phase-amplitude results in Figure 26, is 
then consistent with the previously noted departures from 
potential flow theory. 

D. PERIOD 4 • 

The significant features noted in this last period are 
a small negative u-w correlation, a shift away from near 
potential flow for u and w, and definite stress deviations 
which suggsst wind-wave interactions. 

With respect to momentum transfer processes, r|(u,w) 
patterns in Figures 27 and 28, (Panel C) , appear as minor 
elliptic, asymmetric circles. This appearance signifies 
a smaller probability of joint occurrences of negative 
correlation of u and w. However, the inner rings in 
Figure 27 (Panel C) , represent a slight elliptic pattern 
in the opposite direction and reveal a definite positive 
correlation. Similarly, the phase relations of u and w in 
Figure 29 show u and w to be nearly in-phase. In subset 
3 of this period. Figure 28, (Panel C) , the positive corre- 
lation is not as evident suggesting a shift between u and 
w phases which can be verified by the phase-amplitude re- 
sults in Figure 30. 
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Figure 28. JPDF-CMF Patterns for Period 4, Subset 3. 
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Figure 29. Phase-Amplitude Results (Components) for 
Period A, Subset 1. 
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The orientation of CMF patterns in Figures 27 and 28, 
(Panel B) , corresponding to u (f) , T]) are similar to those 
for near potential flow during the beginning of the period, 
then gradually rotate as the period progresses to a final 
rotation of about 30 degrees behind the crest. Similar 
rotations are not evident in the V7(r], n) CMF patterns, 

(Panel F) . Phase relationships in Figure 29, however, show 
maximum values of w occurring forward of the front node, 
a feature noted by Karaki (1968) and Stewart (1970) . A 
possible reason for this rotation of the CMF pattern may 
be the turbulent interaction between the near surface wind 
and the wave. 

With respect to phase variations of the stress uw, 
initially the maximum positive values of uw , Figure 31 are 
forward of quadrants II and IV and are not equal in magni- 
tude. Later in Figure 33, the magnitude of uw is nearly 
equal in both quadrants II and IV. In both cases, however, 
a sufficient shift is noted from potential flow. This was 
also observed by Kendall (1970). The uw(f|, r\) patterns. 
Figures 27 and 28, (Panel A), also reflects these deviations 
from potential flow since higher CMF values appear in 
quadrants II and IV. However, looking at the phase rela- 
tionships in Figure 31, quadrant IV, downward transfer of 
excess momentum is the largest for both large and total 
waves. During the middle of the period. Figure 32, quad- 
rant II is predominant over large waves, while quadrant IV 
is the^ dominant over total waves. During the last part of 
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Figure 30. 



Phase— Amplitude Results (Components) for 
Period 4, Subset 3. 
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Figure 31. Phase-Amplitude Results (Stresses) for 
Period 4, Subset 1. 
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Figure 32. Phase-Amplitude Results (Stresses) for 
Period 4, Subset 2. 
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Figure 33. Phase-Amplitude Results (Stresses) for 
Period 4, Subset 3. 
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the period, Figure 33, large uw stress amplitudes appear in 
quadrant II for both large and total waves. In general, 
there is a shift from an initial downward transfer of 
excess momentum to an upward transfer of deficit momen- 
tum, This suggests that, due to the waves* influence, 
momentum is being transfered towards the surface during the 
first half of the period and transfered away from the sur- 
face during the last half of the period. During the 
reversal of momentum transfer between the first part of 
the period to the last part of the period, mechanical tur- 
bulence may well be the important factor in the momentum 
transfer . 
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SUMMARY AND CONCLUSIONS 



Analyses were performed on observational data obtained 
over ocean waves. These analyses provided results on the 
phase and amplitude relationship of the wave-induced fluc- 
tuations. These results were interpreted with respect to 
potential flow predictions for a wave propagating beneath a 
shear flow and the role of turbulence in wind-wave coupling, 

A significant feature was the consistent rotation of 
the CMF pattern in the u(f], X]) and w(f), X]) results from the 
orientation which would correspond to potential flow. This 
rotation was noticeable in all periods analyzed to some de- 
gree. Also observed was the absence of the quadrature rela- 
tionship of u and w as predicted by potential flow theory. 

In comparison to potential flow predictions which have 
equal magnitudes of uw phase-amplitude in all quadrants, 
the uw phase-amplitude results revealed variations of magni- 
tude over different phases of the wave. This suggests that 
nonlinear processes are responsible for momentum transfer 
in wind-wave coupling. 

Therefore, wave-caused motion of the air in the near- 
water layer was observed to be essentially nonpotential. 

This was most evident in periods 1, 3, and 4 with a lesser 
deviation from potential flow during period 2. However, 
during period 2, a marked change was taking place in the 
physical make-up of the wave. 
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